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The drag of Ineltalled aircraft radlA'tor's has "been de- 
termined In various caaee In wind-tunnel and free-flight 
tests (references 1, 2, azwl 4). .Unf ortunateljf hoverer, 
the magnitude of the total drag alone affords no informa- 
tion about the quality of the l:'la,dlator Installation* The 
present report therefore proposes the determination of the 
ahsolute magnitude of the' total radiator drag and. In ad- 
dition, of the different causes of the radiator drag. 

In figure 1 is shown a normal radiator installation 
for a liq.uld-cooled engine. Now assume the radiator in- 
cluding cowling removed from the nacelle and mounted into 
a frlctlonless flat plane exposed to an air current as in 
figure 2. The contours In direct proximity of the radiator 
installation shall remain unchanged. The -eventual differ- 
encos in drag and flow volume resulting from the two in- 
stallation arrangements are to he attrlhutahle to the mutual 
interference "between radiator and -airplane components. 
Since this effect is not dealt with here, the analysis can 
be restricted to the simplified case of installation of 
fj.gure 2. 

Two sones, A and B, are to he recognised in the field 
of flow of the radiator: zone A to Include all streamlines 
paseing through the radiator, zone B all remaining stream- 
lines and up to the plane of intersection a-a passing 
through the cooling air discharge orifice (fig. 2). The 
boundary' of the two zones forms an area oompoeed of stream- 
lines which partly coincides with' the contour of the radi- 
ator cowling. Visualizing this surface as a fixed wall, the 
total drag of the. radiator arrangement can be divided into 
drag Vj^ and V3, whereby drag is equal to the force 
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In Inflow direction on.tlie 'boundar7 of zone. A and drag 
tf^ equal to the force, on. the "bound ary of zone B. 

■ ■ dotation 

T, Telocltles In the. plane of Intersection a-a. 
Vg^t velocity in cooling-alx-discharge section, 
v^i ▼oloclty in asstunodly free radiator section. 
' 7q, Telocity in -andlsturhed flov. 

p, .a'bsol'ate pressures in plane a-a. 
Pa'- 'T?'^es\MTo in coollng-alr-diachargo opening. 
Fj^i ' frontal area of radiator element. 
■(f,..tota^ drag of intstallod radiator. 
^a%-- ^^^a of cooilng-air-dlschargo opening. 

p, air ..density ■ 

^q's 7o^^. dy.naffl.i.c pressure of air stream. 

7, 'cooling-air volume in time unit. 
Ap^t pressure gradient in. radiator element. 

'Subscripts X, y» z, denote the components of the • 
velocity in the coordinate directions .z, y, z. 

If assuming the flow devoid of any sourpe or sink, 
the drag iif^ and Hf^ follows from the Impulse theorem 
dt 

,B 



= yy -p + p ▼xC'o— (1) 

. . ■ • A : . ■ 

l^A = J^J^ -P + P Vx(vp -'Vi)..d]P ■ ■ (2) 

if = Vj, + = °w I IPk ■ (3) 
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where' dlT dealgnates "A etlrfeloe element ^ -euld tbe Integra- 
tion ••eat-enda osr^'.-tlie' s.one ■ -A or B 'of plan^ ■ . a-a. When 

e " P + ■§ (▼x'* + ^y^'+ = % ▼o* - As (4) 

where g Is the total pressure and ^g.'the energy Iobb 
of -a Btre'aA filament up to plane ' a-a "relative to a-s'tream 
filament at Infinity V^^ ■ can equally he vrltt'en as 



B 2 



rr [-(vo - + vy* + ▼b^Io + yy^ ^g di (5) 



= 2 ["(^° - Vx)* + ^y' + + yy Ag dF (6) 

Now the assumption Is nade that the pressure and the 
speed In the coollng-alr-dlscharge opening can he expressed 
hy a mean value and furthermore • that the velocity compo- 
nents along the y, ss axis In the qoollng-alr-dlsehaz'go 
opening are negllgihly small. Then the relation for rf^ 
hecomes 

*A'f [^Sa- ft'o-r^j'Jira " 

fig^ ' Indicates the energy loss on passing through the radi- 
ator and amounts to 



■'.A«a Pa - f .V + Ito' 

Thus with' given floW v.olume and known "pros sure In the cool- 
lng-alr-dlscharge opening the'drdg -Vit can he computed. 

As concerns the magnitude of V^, the following pret- 
dlction cen be made; with aerodynamic ally correct radiator 
cowlingi the energy losses in ■ circulation ahout the radia- 
tor vlll he quite loWf so that the second term of equation 
(6) is, in this case, ne-gllglhle. The. magnitude of the 
first term can he estimated as follows Tlsualise the radi- 
ator replaced hy a certain source and sink 4iB-'^'i^'0''l>i oh in 
Bone A, which is so arranged that the boundary surface 
of zones A and B beeoQcs 9 stream surface.' Ifhe force. 
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exerted Internally on the stream 'eurfaoe mustf in this 
caaec "be equal t.o the pressure exerted on. the stream sur- 
face from the outside « I.e., In this case 

Assuming that the speed in the discharge opening is replace 
ai>>a hy the average Ta-lue Ts* Bind, corresponds to the pres 
sure Pa hy los.8~free flov. the yield of the source- 

sink dlstrlhutlon must assume the following yaluet 



(^a* - ''a) 'a 



where 



I- ▼a" =f 'o** - Pa (10) 
The momentum theorem applied to zone A gives 

^B + P ^o <i = P ▼a*^ 'a - P ^a ^o. ^^a + Pa ^a ^^^^ 

Wb = I (^o - va')^ I-a (12) 



(13) 



Incidental to the application of the momentum theory, the 
propulsion of magnitude p "Vq Q exerted on the. sources 
and sinks must, of course, he also taken Into account. 

Admittedly there might he cases where the above as- 
sumptions are no longer permlssihle and other sources and 
sinks outside of sone A oust be assumed. In such cases 
the expression, equation (13) assumes, as a rule, smaller 
or even negative values, but not exceeding the amount ac- 
cording to. equation (13). (Compare, for Instance, the as- 
sumptions acoordlnjg to (3), (5).) Since the pressure p^^ 
nearly always ranges botvecn 0 and ±0.5 qof the expres- 
sion according to equation (13) always remains small rela- 
tive to and can therefore be disregarded. 

With equations (3), (7), (6), and (13) « and given 
dimensions, known cooling volume and known prescfure in the 
coollng-alr-outlet opening, the air resistance of the radi- 
ator assembly can be predicted. This resistance computed 
on the assumption of loss-free flow outside the radiator 
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cowling -iBk'- In contradistinction from the neaaiired tojbal 
drag. V, designated w.ltb. V*. If there' Is a difference 
"between- tf and- V*, 'it mxtB-t he due to energy losbas out- 
side the radiator* I.e., the' second term of equation (5) 
will not -he equal to sero. Therefore: 



W = W* + / / - g dl = c^ | vj,3 Ijf. ■ (14) 



= Ow* f ▼o* Ti •= Wa± • 

In practicet of course, the avoidance of energy losses out- 
side the radiator oowllng, say, hy greater snoothnoss of 
the cowling, etc., wiir'nake the changes in and Pg^ 

quite Boall, so that tho difference of W and tf* does 
not exactly correspond to the actually ohtalna'ble gain in 
drag'. However, this effect Is very small and disregarded 
for reasons of simplicity. 

The radiator drag.. Is chiefly, due to the energy losses 
on flowing through the radiator. Xo con'bat the- power in- 
put, given hy coollng-alr volume T and the pressure gra- 
dient Ap2 in the radiator element, a certain mialnum drag 

m ^ "^o ^ K Wjj, 2 o s 

must he taken in . the hargaln. 

■ 

But in reality, tho energy loss Ag^ at passage is 
much higher t'han corresponds to the pressure gradient Ap. 
in the radiator element. Tlgu^ed with this value. It af- 
fords a drag' Wj^s . 

Wl = V Aga A- = Ag^ Ij. |£ = c^j I Vo= (17) 



The differences hetween Vq and Vj are due to flow 
losses on passage through the radiator cowling, the differ- 
ences "between and the computed W* due to the hasic 
disposition of the radiator* The detomlnation of V*, 
Vq* tf].! and V afford an insight into the composition of 
the radiator drag. 



6 
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V, V*, Vq, Vj. can e(iiLall7 "be appraised mathemat- 
ically If 'based on the folldvlng simplifying aesxiaptlonBa 
Limited to a two-dimensional prohlem the flow past the 
radiator Is Tlsuallzed as' 'being 'bounded 'by a flat wall, 
as In figure 4, arranged, parallel to the plane of Instal- 
lation. The else of the- tunnel section V Is put In a 
certain ratio to the radiator dimensions and so determined 
that the oondltlonB at the coollng-alr outlet agree as 
much as possl'ble with actual conditions. Since section J 
does not exceed a certain slzsi the Telocity outside of 
the coollng-alr-outlet section In the plane' a-a can 'be 
equally replaced by an arerage value t^. The pressure 
In this section Is assumed constEwt and must be equal to 
Pel. From the equation of continuity follows: 

J'a " 

(vi - Vo)(3P - IP J = (Vq - Tg) (18) 
P F 

Y. = v_ - — — — — T™- 

^ P-Pa P-Pa 

Assuming zero pressure in the plane e-o directly be- 
fore the radiator, and designating with Agg the energy 

'losses at passage, with ts± the energy losses at flow 
around the radiator, from plane e-e to plane a-a, we 
have, according to equation (4): 



Pa = 4o - f ^a"* - Aga (19) 
Pa = 4o - 2 ^1* - ^61 ^20) 

Once the flow losses Lg^ euid Aga have been ascertained, 
the cooling volume Euid the pressure In the outlet opening 
can be computed with eqtiatlons (18), (19), and (30). Aga 
consists of the pressure loss Ap^ In the actual radiator 

element and the loss at Inflow and outflow. 

For Ag^ and Agj^, the following appropriate assump- 
tions are made: 

Aei = ^1 f ^o' (21) 

Aga = Apj + £a I ^o^ ^22) 
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^^Z." a ""o' ^^^^ 

■ - ^ =■ " ... . ■ .. , . ■ 

^li t&* ^^'^ "^^^ conatanta. ' She" adnlSBi'billty of 

the afmplif icatlona effected through equatlona '(21-)« (22), 
and.. (23) . miiat tb ohecked In th9 Indlvldua}. oaaea, or elae 
properly modified reiatlobs Introduced for' Ag^. Ag'^, 

and ' Ap^* In the following caloulatlona , the noAdlmenslonal 

brr Pa ^StL 

quantities, — . -■ , and the nondlaenslonal ooef- 

. 7o lo lo 

flcienta, Cy*, Oy, Cy^, Cy^; oonf orma'bly to equations (3), 
(15), (16), and (1?) are employed. 

The application of equations (15) and (17) to the 
present case gives: 

W = /Sga3'a+Agl(r-I'a) -■|(vo-Ta)^S'a-f (▼o-^l)^(T-I'a) (24) 

■■ TT-n^ ' ^'^^) (.-pt; ^^^^ 



* *o i^K \ \ . ■'o/ V r / 

OE ^Sa 



^ ^0 lo 



(27) 



cw (28) 

Ag_ v_ tj Ag4 
with the quantities, , — , — , — to 'be detarmlned 

<lo '^o <lo 

according to equations (19) to. (23).. 

Now the derlvated equ'atl'ons are applied to an Illus- 
trative example. The selection of the constants and the 
geometric' dimensions Is as follows: 

■ ■ ' ' >K ■ ' 

■< — = 0.3'2,' Cw_ = 6 

tl = .la 



I 
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A flap fitted at the coollng-alr . out let (fig* 4) af- 
fords different ooollng-a'lr-outlet sections T^^. At first, 

the flov Tolume ^ f.or each flap setting Is determined; 

■ ■ "^o ■ ■ • p ■ 

thle Is best done grapjxlaal^y, the value of 1 - — = , for . 

err" *lo- 

each flap sotting "being plotted against — = and according 

■ Vo • ■ 

to .equations (19). and (20). The Intersection of the two 
curves gives the desired flow volume ve/vq* 

In this manner, the flow volume ^ of figure 4 was 

T ■ ■■ ° ■ 
determined for flap setting _ q, 0.25, 0.50, 0.75, and 

the related values , according to equa- 

te 1o • ■'■o 4o 
tlona (18), (19), and (20) aacnrtained, and then the values 
of Cy, c»*, Cy and Cy. computed according to equations 

(35) to (28) and connected "bj curves. She effect of any 
change.B, as, say In the geometric dimensions, the assump- 
•tlons regarding, flow losses, etc., on drag, flov volume 
and pressure In the coollng-alr outlet are readily detected 
hy this method of calculation. 

In the following, .the arguments are applied to several 
experlnental radiator investigations. The radiator instal- 
lation of figure 6 with nacelle cuid a piece of the wing was 
explored at full scale in the wind tunnel* Dimensions and 
arrangemont of radiator installation aro in the principal 
points those used in the other example* She radiator as- 
sem'bly was studied at different flap settings and drag, 
flow volumo, and cooling-air -outlet pressure ascertained* 
■The radiator drag followed as difference from the recorded 
drag less the drag hy insteJLled radiator and host naeolle 
covering. The results are tabulated in table I. 

The drag coefficients c^ , c^ , c^* were determined 

"m " i " 

Vlth the help of the foregoing argument a . and plotted in fig- 
ures 7 and 7a. tfith arrangement 1, a large portion of the 
drag was fouiid to be caused by the flow losses while circu- 
lating around the radiator. Then the radiator inlet vaa 
shortened and the inlet edge rounded off, which gave arrange- 
ment 2, on which these losses, especially at small flap 
settings, could bo materially lowered, although it did not 
affect the losses at passage through the radiator, according 
to figure 7* An Intermediate baffle fitted at the radiator 
inlet conformahly to arrangement 3, figure 7a afforded con- 
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slderable Inprovooout lii.'.'.."Oy^ - au-d hsnce Indirectly In the 

'total drag. Ihle exanple Indl&ates how. vlth.the aid of 
this aethodf the Individual oaupes of the radiator drag can 
ho ascertained. A' s.li!illar iStiTeBtlg&tldn at smaller scale 
had heen nade earlier' for a different radiator' lAstallatlon. 

A coaparieon af the c caput .e.d curves of figure 5 and 
the measured oury.eiB of figure 7 Indicates' Very cloee agroe- 
BLont, irhloh,. oonalder.lng tho extent of neoesBary slnpllfi- 
catlons 1b Burprlalng. The a.B'Bunptl on. regarding the ex- 
ternal flov losaes ls"haturall7 o'niy d'osditTonally correct i 
since* with- flap closed, these losses .nust unAouhtedly he 
Biihstant ially greater than with 'f!id,p. open, 'a fact suhstan- 
tliated by .experiment. The sane also holds for flow losses 
through the radiator 'cowling. 

By. neans. of tho foregoing arguncnts the prohlea of 
radiator drag reduces to the solution. of a series of par- 
tial prohleas each of w'hich can he dealt with separately. 
It affords a clear picture of the coaposltlon of the radia- 
tor drag and a predoternlnatl on ' of tho radiator dreg and 
cooling capacity In most' cases with sufficient accuracy. 

Tho enployed method is equally applicable to air- 
cooled engines. 



Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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• ' ' Table -I -'- 

Test Heaulte vlth Sadidtor Arraugement of ]rigur« 4 



'a 


▼o 


2a 


Pa 




ao 




(ISO) 






0.000 












0.320 








.148 


0.104 


0.700 


-0.06 


0.57 


0.119 


.234 


0.071 


0.059 


0.0123 


1 .281 


.182 


.650 


-.04 


.618 


.222 


.234 


.139 


.113 


.0403 


.529 


.306 


.580 


-.235 


.899 


.537 


.401 


.382 


.275 


.1643 


.750 

• 


.355 


.473 


-.383 


1.161 


.68? 


.740 


.662 


.412 


.244 


0.000 












0.162 








.148 


0.100 


0.675 


+0.01 


0.555 


0.113 


.106 


0.066 


0.055 


0.0113 


2 .281 


.179 


.636 


-.03 


.626 


.214 


.145 


.139 


.112 


.0383 


.515 


.314 


.608 


-.24 


.871 


.563 


.380 


.369 


.874 


.1767 


.750 


.390 


.520 


-.41 


1.139 


.815 


.662 


.681 


.444 


.318 


0.000 












0.179 








.148 


0.104 










.086 








3 .284 


.212 


0.747 


-0.03 


0.475 


0.280 


.126 


0.116 


0.100 


0.0594 


.750 


.370 


.493 


-.395 


1.152 


,739 


.665 


.672 


,426 


.2734 


.529 


.306 










.360 
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figure l.~ General arrangement of a 
radiator aasembly: for a 
liquid-cooled engine, mounted in the 
engine micelle of an airplane. 
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Figure 2.- Radiator installation in infinitely extended frictionlese 
plane. 
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along plane a-a 

Zone f) 

////// 
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Figure S^- A specific source and sink distribution substituting for 
the radiator. 
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Figure 4.- Radiator in a closed chaimel 




Figs. 5,6,7,7a 
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Figure 5.- Theoretical de-termination of radiator drag and flow volume. 

Direction of F/fght 




Arrangement / and S tvithout intermedia^a 

Arrangement 3 corraaponds to arrangement S 
tri'fh intermediate baffle 

Figure 6.- Radiator mounted on engine nacelle. 

Arrtmgement S 

r : <7c 07 

Arrangement I 
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Figures 7 and 7a;- Experimental radiator drag and flow volume. 



